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ABSTRACT 

Four techniques for detecting multilayered clouds and retrieving the cloud properties using satellite data are explored to 
help address the need for better quantification of cloud vertical structure. A new technique was developed using 
multispectral imager data with secondary imager products (infrared brightness temperature differences, BTD). The other 
methods examined here use atmospheric sounding data (CCUslicing, C02), BTD, or microwave data. The C02 and 
BTD methods are limited to optically thin cirrus over low clouds, while the MWR methods are limited to ocean areas 
only. This paper explores the use of the BTD and C02 methods as applied to Moderate Resolution Imaging 
Spectroradiometer (MODIS) and Advanced Microwave Scanning Radiometer EOS (AMSR-E) data taken from the 
Aqua satellite over ocean surfaces. Cloud properties derived from MODIS data for the Clouds and the Earth’s Radiant 
Energy System (CERES) Project are used to classify cloud phase and optical properties. The preliminary results focus 
on a MODIS image taken off the Uruguayan coast. The combined MW visible infrared (MVI) method is assumed to be 
the reference for detecting multilayered ice-over-water clouds. The BTD and C02 techniques accurately match the MVI 
classifications in only 51 and 41% of the cases, respectively. Much additional study is need to determine the 
uncertainties in the MVI method and to analyze many more overlapped cloud scenes. 
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1. INTRODUCTION 

Improved characterization of the vertical structure of clouds is needed for the development of better cloud 
parameterizations and the validation of cloud fields in climate models. Satellite data constitute the most important 
source for quantitatively describing the distribution of clouds at large and small scales over extended areas. 
Traditionally, analyses of satellite data have been limited to interpreting the radiances in terms of single-layer clouds 
resulting in large errors in the retrieved properties for multilayered cloud systems. Recently, a variety of methods have 
been developed to detect multilayered clouds and to estimate the properties of those clouds classified as multi-layered. 
Combining visible (VIS; 0.65 pm) and infrared (IR) retrievals of cloud properties with microwave (MW) retrievals of 
cloud water temperature Tw and liquid water path LWP appears to be a fruitful approach for detecting and retrieving 
overlapped clouds 1 ' 2 ' 3 . A method using 1.6-pm reflectance and 11-pm brightness temperatures, T, shows some potential 
for unscrambling the cloud impacts on observed radiances 4 , but it has not been extensively used. The use of BTD, the 
11-12 pm brightness temperature difference, in conjunction with x, the retrieved visible optical depth, was suggested 5 as 
a means to detect multilayered clouds. A slightly different approach was developed using the BTD along with the VIS 
reflectance 6 and has been applied to detect multilayered clouds using daytime data taken over most non-polar regions 7 . 
The discrepancy between cloud-top pressure derived from a CCE-slicing retrieval and the IR-based cloud pressure has 
been exploited to detect overlapped clouds and retrieve the properties of each layer 8 over a large portion of the Earth 9 . 
These recent ventures into passive remote sensing of multilayered clouds are encouraging, but the accuracy of these 
retrievals and their limitations, especially regarding detectability, are poorly understood. To begin characterizing which 
clouds are detected and retrieved properly, this paper performs some initial comparisons of these techniques and another 
method recently developed for application to both geostationary and polar-orbiting satellite imagery. 
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2. DATA 


Data from the Moderate Resolution Imaging Spectroradiometer (MODIS) and the Advanced Scanning Microwave 
Radiometer-EOS (AMSR-E) on the Aqua satellite are used to perform the various retrievals. The 1-km MODIS data are 
analyzed for the Clouds and the Earth’s Radiant Energy System 10 (CERES) using the VIS-IR-Solar-infrared-Split- 
window Technique" (VISST) to retrieve single-layer (SL) cloud properties for each pixel 12 . The CERES cloud retrieval 
output, which includes the cloud properties such as cloud optical depth x and cloud effective temperature Tc, as well as 
the original pixel radiances at 0.64, 1.6, 2.1, 3.7, 11, and 12 pm, for each pixel is combined with the radiances for the 
1.3, 4.0, 6.7, 8.5, 13.3, 13.5, 13.9, and 14.2-pm MODIS channels. Cloud LWP and Tw are derived from the AMSR-E 
MW data at the 12-km resolution 36.5-GElz field of view (FOV) by matching the multispectral MW data to radiative 
transfer model (RTM) calculations 13 . The CERES pixel-level results are then convolved into the AMSR-E footprints 
following an earlier approach 2 in order to apply the MW-VIS-IR (MVI) multilayered cloud detection method 1 and 
retrieve cloud properties using the VIS two-layer cloud retrieval technique’ VTRT. In this pilot study, the CERES 
MODIS pixels are matched to the nominal nadir-size footprint of the AMSR-E to facilitate comparison of the different 
techniques. 

The comparative analyses focus on one granule of Aqua MODIS data taken 13 July 2004 off the coast of Uruguay at 
1316 UTC. The RGB MODIS image in Figure 1 shows a complex system with a mixture of ice (pinkish, grayish), 
liquid water (white, yellow, and peach), and overlapped clouds. Also shown are the VISS-retrieved values of x for all 
cloud types assuming SL clouds. The cloud optical depths range from 0.3 up to 128. 



Figure 1: CERES imagery for multilayer cloud detection comparisons for 1316 UTC, 13 July 2004. left: Aqua MODIS pseudo-color 
RGB image, right: VISST-derived cloud optical depth. 


3. METHODOLOGIES 

Four different ML cloud techniques are compared here; the MVI, the VIS BTD method 6 (VBM), the CCVslicing 
method 8 (COM), and the CERES BTD method (CBM). The MVI is used as the reference technique because it can 
detect water clouds over ocean in the presence of ice clouds as long as no precipitation occurs. While it is not possible 
to exactly duplicate the published techniques, every effort is made here to match them as closely as possible. 

3.1 MVI 

The MVI is used only if the VISST retrieves at least 98% ice cloud coverage for the MODIS pixels in the convolved 
AMSR field of view and the solar zenith angle is less than 78°. The scene is classified as being composed of SL ice 
clouds, if the retrieved MW LWP < 40 gm" 2 . If LWP > 40 gm" 2 , Tw - Tc > 5 K, and no precipitation is detected, the 
scene is classified as ML ice clouds, that is, an ice over a water cloud system. The precipitation detection is based on the 



amount of polarization in the 36.5-GHz channel. Single-layer water, clear, partly cloudy, and mixed phase scenes are 
not considered here. 

3.2 VBM 

The VBM has two formulations, one for an imager with only a few channels and one with many spectral channels like 
the MODIS 6 . The multispectral approach, approximated here, uses BTD, the VIS reflectance R v , T, and the 1.3 and 1.6- 

pm reflectances, R , and R n , respectively. The pixel is first assumed to be SL. It is reclassified as ML if either one of the 

following tests are satisfied over ocean, snow, or ice surfaces. 

3.2.1) 0.3 < R v < 1.0 and BTD > THb and T < 270 K. 

3 .2.11) 0.025 < Ri < 0.4 and R n > THn and T < 280 K and R„ / R r < 1.0 and (BTD < THb , if R, < 0.08). 

Values for the BTD and near-infrared thresholds, THb and THn, respectively, are available in lookup tables 6 . THb 
depends on surface type, R v and the viewing and solar zenith angles, while THn varies as a polynomial function of /?,. 
The Aqua MODIS 1.6-pm channel has not functioned properly and is essentially useless. Since it is similar in 
reflectance behavior over ocean and clouds, the 2.1-pm channel is used here instead of 1.6 pm. To account for the 
differences, the 1.6 and 2.1 -pm channels on the Terra MODIS were compared over clouds to derive a near-infrared 
threshold THn(2. 1 ) in VBM test 3. 2. II. This initial threshold value is 

THn(2A) = 0.58 THn(\.6). (1) 

The validity of this threshold formula can only be determined by further application of the VBM to Terra MODIS data 
using the 1.6 and 2.1-pm channels separately for the same pixels. 

3.3 COM 

The COM method is applied in the following manner with the initial assumption that the pixel contains only SL clouds. 
The following tests are applied to pixels having a valid VISST retrieval as a an ice or water cloud. For water clouds, the 
pixel contains multilayered cloud if 

3.3.1) e(V) - e(C) > 0.3 and z( C ) - z(V) > 1 .5 km and e(V) > 0.9 and r e > 9 pm and p(C) < 450 hPa, 

In this test, e, z, and p denote the cloud emissivity, height, and pressure, respectively, as computed from the VISST (V) 
or the C02-slicing (C) methods. The VISST-derived effective droplet radius is r e . 

For ice clouds, the optical depth limit is computed as 

ODL = 0.96 -1.09 [pln[l -s(C)], (2) 

where the cosine of the viewing zenith angle is p. The cloud is classified as ML ice is the following test is satisfied. 

3.3 .11) p( C) < 500 hPa and x(V) > ODL and e(C) < 0.85. 

An additional set of tests is applied if the cloud is ice and x(V) > 20. The cloud layering is classified as indeterminate if 
any one of the following statements is true. 

3. 3. Ilia) -0.5 K < BTD < 0.5 K. 

3.3.IIIb)3.0K< T- T( 3.7 pm) <3.0 K. 

3.3.111c) -3.0 K < T- 7(4.0 pm) < 3.0 K. 

3.3.IIId) -3.0 K <7 - 7(6.7 pm) < 3.0 K. 

3.3.IIIe) -0.5 K< 7 - 7(8.5 pm) < 0.5 K. 

3.3.IIII) -3.0 K< 7- 7(13.3 pm) < 3.0 K. 



The values of the brightness temperatures T(k) correspond to observed value for the given channel wavelength, A. The 
tests in 3. 3. Ill are designed to account for the uncertainty in the height of the lower cloud and are not part of the original 
COM 9 . Whenever a thin high cloud occurs above a thick ice cloud, it can produce brightness temperature differences 
approximately within the ranges given in 3. 3. Ill limits. The indeterminate category basically assumes that the ice cloud 
is too thick to reliably determine whether a water cloud resides underneath the ice cloud. This type of condition often 
arises for convective clouds. 

The COM serves as a starting point for a retrieval process that first searches for the average properties of nearby low 
clouds. They serve as the initial properties of the lower cloud, indicated by the subscript /, for use in RTM calculations 
of the VIS reflectance at the satellite view. The cloud input for the RTM consists of %\, r e \, and 7) from the neighboring 
pixels for the lower cloud and x u , D eu , and T u from the C0 2 -slicing and VISST retrieval for the upper cloud. The C0 2 - 
slicing retrieval and lower cloud properties are iterated with adjustment of the appropriate parameters before each new 
step until the calcualted VIS reflectance matches R v . The output includes cloud properties for each layer. 

3.4 CBM 

The CBM is a technique, newly formulated following the suggestion of earlier research 5 , that uses the VISST-derived 
parameters and the BTD to classify cloudy pixels as SL, ML, or indeterminate. Figure 2 shows the flow diagram for the 
CBM. It begins with the VISST retrieval and applies different tests depending on the phase of the SL retrieval. In 
addition to x, the CBM uses either the effective ice crystal diameter D e or r e depending on the phase. In addition to the 
tests listed in the flowchart, the CBM applies the same indeterminate tests (3. 3. Ill) as the COM for all ice clouds with x 
> 20. For clouds classified as ML, the upper and lower cloud properties are estimated in manner similar to the COM, 
except that the Solar-infrared Infrared Split-window Technique" (SIST) is used to estimate x u because no C0 2 bands 
are used. To improve the height determination, the SIST is also used for SL clouds when the VISST retrieves x < 4. No 
iteration is performed in this prototype version. 
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Figure 2: Flow chart for classifying VISST retrievals as single or multi-layered clouds in the CERES BTD Method. 


a) Cloud Physics Lidar, ER2 / CRYSTALFACE, July 26, 2002 
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c) CERES Multilayer Algorithm Results 
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Figure 4: Multilayered cloud case over the Gulf of Mexico with matched MODIS and ER-2 CPL data, 1830 UTC, 26 July 2002. 


An example of the CBM applied to Aqua MODIS data is shown in Figure 3 for a match with the Cloud Physics Lidar 14 
(CPL) on the ER-2 as it flew over the Gulf of Mexico during 26 July 2002. The CPL detects a thin cirrus layer (Figure 
3a) between 13 and 15 km above scattered cumulus clouds between 0.5 and 2 km. The BTD (Figure 3b) along the flight 
path varies between 2 and 4 K with larger values coinciding with slightly thicker cirrus and the absence of cumulus 
clouds. The CBM classification (Figure 3c) detects ML clouds only when the lower cloud appears to be thick (bright) in 
the CPL return. Although some ML clouds are misclassified as SL, most of the SL clouds are accurately designated. 
The retrieval method yields mixed results for cloud heights (Figure 3d). At 18.3 UTC, it only produces a slight 
separation between the high (red) and low (blue) cloud, possibly because the low cloud height was incorrectly specified 
as being too high. At 18.37 UTC, the clouds are correctly classified as SL and the VISST (black) places the cloud 
between 9 and 10 km. The effective cloud height from VISST should be located somewhere between the cloud base and 
top. The SIST yields a more accurate retrieval of Z, placing the cloud between 13 and 14 km. At 18.40 UTC, the clouds 
are categorized as ML. The VISST retrieval is basically the same as the low cloud (blue), while the high cloud retrieval 
is between 9 and 14 km, certainly an improvement over the VISST retrieval. At 18.42 UTC, the ML cloudiness is 
misclassified as SL because the lower cloud is optically thin. In this case, the VISST retrieval is more accurate than the 
SIST. Results for the remainder of track are similar to those from the first part of the flight. This prototype retrieval is 
promising, but the iterations and conditions for application need to be implemented to fully assess its capabilities. 




Figure 4: Cloud layering classifications for scene in Figure 1. (a) MVI: blue - precipitation, yellow - ML ice, magenta - SL ice. (b) 
VBM: white - SL ice, gray - SL water, yellow- ML ice from 2.2.1, red - ML ice from 3.2.11, dark blue - ML water from 3.2.1, light 
blue - ML water from 3. 2. II, green - clear. 


4. RESULTS AND DISCUSSION 

The four ML detection methods were applied to the scene in Figure 1 and are compared below. Figure 4 shows the 
classifications from the MVI (Figure 4a) and the VBM (Figure 4b). The MVI is restricted to the center of the swath 
because the AMSR-E scans conically at a viewing zenith angle of 55°. Thus, all comparisons with the MVI use only 
those MODIS pixels over water taken at a viewing zenith angle of 55° or less. Results from the COM and CBM are 
shown in Figures 5a and 5b, respectively. The MVI precipitating clouds correspond to indeterminate, ML, and SL ice 
clouds from the other methods as well as to some SL water clouds. The MVI ML clouds typically correspond to 
different ML clouds from the various other techniques. The CBM and VBM are most similar, at least, in the yellow 
areas, however, the VM 3. 2. II test (red) classifies more pixels as ML. Some correspond to MVI ML pixels, while others 
coincide with MVI SL pixels. 



Figure 5: Cloud layering classifications for scene in Figure 1. (a) COM: brown - indeterminate, yellow - ML ice, white - SL ice. (b) 
CBM: brown - indeterminate, white - SL ice, gray - SL water, yellow- ML ice, dark blue - ML water, green - clear. 



Table 1. BDM classification of MODIS pixels in marine AMSR-E pixels containing 100% cloud cover as determined by the VISST 


for the same MODIS pixels. Results in % of total number of MODIS pixels, except N(AMSR). 


MVI class 

SL water 

SL Ice 

ML water 

ML ice 

indeterminate 

clear 

N(AMSR) 

ML ice 

0.4 

27.2 

0.1 

29.3 

n/a 

0.0 

11632 

SL ice 

0.1 

9.9 

0.0 

11.2 

n/a 

0.1 

4429 

precipitation 

2.4 

9.6 

0.6 

8.8 

n/a 

0.1 

4260 

Table 2. Same as Table 1, except for the COM results. 

MVI class 

SL water 

SL Ice 

ML water 

ML ice 

indeterminate 

no high 
cloud 

N(AMSR) 

ML ice 

0.2 

28.2 

0.0 

18.7 

2.4 

7.5 

11632 

SL ice 

0.1 

13.0 

0.0 

6.0 

0.1 

2.2 

4429 

precipitation 

0.5 

5.8 

0.0 

3.1 

7.1 

5.0 

4260 

Table 3. Same as Table 1, except for CBM results. 

MVI class 

SL water 

SL Ice 

ML water 

ML ice 

indeterminate 

clear 

N(AMSR) 

ML ice 

0.3 

30.3 

0.2 

24.2 

2.0 

0.0 

11632 

SL ice 

0.1 

15.0 

0.1 

5.9 

0.2 

0.1 

4429 

precipitation 

2.7 

4.7 

0.3 

6.7 

7.0 

0.1 

4260 


The comparisons with the MVI results are quantified in Tables 1-3. The numbers represent the percent of all MODIS 
pixels within the FOVs of all AMSR-E pixels having an MVI classification. Because the MVI is only used for AMSR-E 
pixels having at least 98% ice cloud coverage from the VISST or are precipitating, the percentages will exclude most 
clear and SL and ML water pixels from consideration. The VBM (Table 1) yields the closest match for ML ice at 
29.3%. It is followed by the CBM (Table 3) and then the COM (Table 2). The latter finds no high clouds in 13% of the 
MVI ML cases. The number of indeterminate pixels is similar for both the COM and CBM because nearly identical 
methods are used for that category. The CBM yields the best agreement with the MVI for SL ice at 15% followed by the 
COM and VBM, respectively. The retrievals for MVI precipitating clouds (PRCP) vary considerably for each of 
methods. The CBM and VBM classify roughly 12% of the PRCP as SL water clouds because both rely on the VISST to 
determine phase. The indeterminate category is assigned to 33% of the PRCP pixels by both the COM and CBM . 
Because it has no indeterminate class, the VBM assigns about 30% of the indeterminate pixels to ML ice and the 
remainder to SL ice. 

An overall score for these techniques relative to the MVI retrievals can be assessed by considering only the non- 
precipitating AMSR-E pixels. A correct classification is assumed if the technique yields ML clouds when the MVI 
detects ML ice and SL clouds when the MVI retrieves SL water. Given that basis, the VBM, COM, and CBM yield 
accurate classifications in 50.3, 40.5, and 50.6% of the cases. The accuracies for the VBM and CBM are nearly the 
same despite the differences in ML classifications. Accepting the assumption that the MVI is accurate, these results 
suggest that the three other ML detection methods are not demonstrating the kind of skill required for reliable detection 
of ML clouds. Some of the errors in the ML methods may be due to inaccuracies in matching the MODIS pixels used in 
the MODIS-only retrievals because the nadir FOV size was assumed for the matching. However, a comparison of the 
images in Figures 4 and 5 indicate that this not likely to explain more than a few percent of the differences in 
classifications. Another source of error is the assumption that in an AMSR-E FOV, the entire area is covered by a water 
cloud when LWP > 40 grn' 2 and by no water when LWP < 40 grn' 2 . This is a reasonable assumption given that the 
detectability of LWP from MW data is highly uncertain when the LWP is below that limit and it corresponds to a lower 
cloud optical depth of~5 for r e = 12 pm. Detection of ML clouds with the MODIS data alone is likely to be difficult if 
x u < 5. However, if the clouds under the are scattered and have optical depths exceeding 5, then the MVI might classify 



the cloud as SL ice while the MODIS methods could classify a third or more of the pixels as ML and be correct. 
Conversely, the AMSR-E pixel could be half covered by optically thick cumulus clouds and be classified as ML ice 
while only half of the MODIS pixels are categorized as ML and be correct. The limitations of using the MVI clearly 
need further examination. 


5. CONCLUDING REMARKS 

Preliminary comparisons of four different multi-layered cloud detection methods confirm that passive remote sensing of 
multilayered cloud systems is a difficult and possibly intractable approach for a variety of cloud systems. No pessimistic 
conclusions should be drawn from the results, however, because only one scene was analyzed, the various techniques 
have not been optimized, and the uncertainties and limitations of the reference MVI method have not been clearly 
defined. A framework for performing more detailed and robust comparisons has been established, thus, facilitating the 
analysis of multilayered clouds over the global oceans and for incorporating the optimized algorithms for both 
published and unpublished detection methods. By further detailed examination of the MVI in conjunction with surface, 
airborne, and satellite lidar-radar systems, it should be possible to determine the uncertainties in a given retrieval. The 
capabilities of the other methods can then be assessed using both the active sensors and the MVI to cover more cloud 
types and viewing angles than would be possible using the active sensors alone. 
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